Collective excitations and instabilities in multi-layer stacks of dipolar condensates 
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We analyze theoretically the collective mode dispersion in multi-layer stacks of two dimensional 
dipolar condensates and find a strong enhancement of the roton instability. We discuss the interplay 
between the dynamical instability and roton softening for moving condensates. We use our results 
to analyze the decoherence rate of Bloch oscillations for systems in which the s-wave scattering 
length is tuned close to zero using Feshbach resonance. Our results are in qualitative agreement 
with recent experiments of Fattori et al. on "^^K atoms. 



Recently quantum degenerate gases and Bose con- 
densates with dipolar interactions attracted considerable 
attention both theoretically [1,2,3] and experimentally 
lii. The long range character and anisotropy of dipo- 
lar interactions is expected to lead to a variety of exotic 
many-body ground states [1,2,7,8] and unusual collective 
excitations [3(. A recent addition to this class of sys- 
tems are ultracold atoms in which the s-wave scattering 
length can be tuned to zero using Feshbach resonance 
and thus enhancing the role of dipolar interactions [4] . In 
particular, groups in Florence and Innsbruck used one- 
dimensional (ID) optical lattices to create multi-layer 
stacks of pancake condensates in the regime of small 5- 
wave scattering length and experimentally studied Bloch 
oscillations (BO) [9IIiq|. In this paper we analyze the- 
oretically collective excitations in dipolar superfluids in 
the presence of a ID periodic potential (Fig. [TJa)). Com- 
pared to the single layer case we find that the roton insta- 
bility is strongly enhanced and occurs at a much shorter 
wavevector, which is determined primarily by the inter- 
layer distance. We also observe interesting interplay be- 
tween roton softening and dynamical instability for mov- 
ing condensates. This implies that both the dynamical 
and roton instabilities should play an important role in 
the dynamics of Bloch oscillations. We compare our re- 
sults for the decoherence rate of BO to experimental re- 
sults of Fattori et al. [9] and find qualitative agreement. 

We consider a multi-layer stack of two-dimensional 
(2D) condensates loaded in a ID optical lattice with the 
magnetic dipolar moments {jUd) oriented either perpen- 
dicular or parallel to the layer planes, depending on the 
direction of the external magnetic field (see Fig. [U (a)). 
We start our analysis with the assumption that layers are 
infinitely large, but will consider finite size effects when 
comparing with the experimental results for the decoher- 
ence rate of BO. We further assume that the ID periodic 
potential (which we take to be in the z-direction) is suffi- 
ciently strong so only the lowest single particle subband 
(wavefunctions, '0j(^), where j = 1,2, •••L is the layer 
index and L is the total number of layers) is occupied 
in each layer, and ^pj{z) can be well-approximated by 
Gaussian functions of width W Q. Integrating out the 
confinement wavefunction in the z direction, we obtain 
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FIG. 1: (a) Dipolar condensate in a ID optical lattice . Ex- 
ternal magnetic field, B, is in the x — z plane with an angle 
with respect to the lattice direction (z). a is the angle 
between the in-plane momentum q of excitations and the x- 
axis. Figures (b) and (c) show collective mode dispersions 
for magnetic fields parallel and perpendicular to the layers 
respectively. The dashed curves are the single layer results 
for a = 0, and the solid curves are for L — 10, 20, and 100 
respectively from bottom to top in (b) and from top to bot- 
tom in (c). The three thin solid lines in (b) are results for 
a = 0, while the thick solid line is for a — n/2 and inde- 
pendent of L. Here we choose as = — O.lSao. (d) Same as 
(c) but with L = 10 and different values of as', as — —0.19, 
-0.38, and — 0.472ao from top to bottom. The dashed line is 
the imaginary part of collective mode for as = — 0.472ao. In- 
set: Momentum distribution of thermal quasi-particles of the 
same system for as = —0.19, -0.38 and -0.47 ao from bottom 
to top with T = 100 nK. 



the intra-layer interaction: 
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with a. and m be- 



ing the s-wave scattering length and the atom mass. 
x{0, a) = cos^ — cos^ a sin^ with being the angle 
between dipole moment and the z-axis (direction of op- 
tical lattice) and a being the angle between in-plane 
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momentum q and the x-axis, which is set to be the 
dipole moment direction when the magnetic field is par- 
allel to the layer (see Fig. [11(a)). Finally, P2{x) = 
|(3x^ — 1) is the Legendre polynomial and F{q) = 
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with Erf(x) being the 



Error function. Similarly, the inter-layer interaction ma- 



trix element between layers j and / (/ = \j 
be also calculated: 
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We note that for the case of perpendicular field (6> = 0), 
V;(q) = -5MMe-|ql^^ as VI^ <C showing a strong at- 
traction at a characteristic length scale. Id. As we discuss 
below, such feature will lead to the wavevector of roton 
excitations being much smaller than the corresponding 
wavevector in a single layer system [3]. 

Our starting point is the microscopic Hamilto- 
nian: H = (2 J(l - cos{kd)) + - /i) 6^ ,6^^, + 

^E/c/ (I^P/c(q)p-/c(-q)Vfe(q), where = q2/2m, 
J is the inter-layer single particle tunneling rate, bc^^k 
is the bosonic field operator for magnetic atoms, and 
Pfe(q) = Ep,/c' ^p+q,/c'+fe^p,/c' is the density operator. 
V{q^^k) = ^^=iVj{q) e~'^^^^ . The chemical potential /i 
is determined from the condition that collective modes 
are gapless in the long wavelength limit. The summation 
on k is within the first Brilluoun zone, —ir/d < k < ir/d. 

When discussing collective excitations we consider a 
general case of dipolar condensate moving in the z- 
direction. Assuming that the moving condensate has a 
macroscopic number of particles in a state with lattice 
momentum /cq, the dispersion of Bogoliubov modes in 
the co-moving frame is given by 

u;{ci,k;ko) = 7{K(fc) + eO(q)] 



(fc)+e°(q)+2nol/(q,fc)]},(3) 



Here E^^{k) = 2Jcos(fcod)(l - cos{kd)), no is the 2D 
condensate density in each layer, k and q are the z- and 
in-plane components of the excitation momentum. In- 
stabilities manifest themselves as imaginary frequencies 
in Eq. (|3]). There are three distinct mechanisms of insta- 
bilities: (i) when kod > 7r/2 and V(0,0) > we find the 
usual dynamical instability [11,12], which can be under- 
stood as coming from the negative effective mass caused 
by the lattice band structure; (ii) when F(0, 0) < we 
have a 3D collapse of the condensate [4]; (iii) instability 
at a finite momentum q due to the softening of roton 
excitations. The latter is a unique feature of long-ranged 
dipolar interactions in low-dimensional systems In 
the rest of this paper, we will analyze distinct roles played 
by these three instabilities in the dynamical properties of 
moving condensates. For numerical calculation, we take 



system parameters similar to Florence's group [9^ on ^^K, 
where the distance between the layers is d = 0.516 jam 
and the strength of the optical potential is = QEr 
(here Er = 4.8/c Hz is the recoil energy). As a result, 
J = 0.0645£^i?, and the Gaussian width of each layer is 
W ~ 0.203(i. Throughout this paper, we will use L = 10 
and no = 10^^ cm~^, unless specified differently. 

To demonstrate the importance of inter-layer coupling, 
in Fig. [TJb) and (c) we show the calculated in-plane col- 
lective mode dispersion of a static {ko = 0) condensate 
for different numbers of layers, L (we set the z-component 
of the excitation momentum to zero, k = 0). We ana- 
lyze cases of the magnetic field (dipole moments) being 
either perpendicular {6 = 0) or parallel (^ = 7r/2) to the 
planes of the condensate pancakes. For the perpendicu- 
lar field case (Fig. [Ifc)), the long wavelength behavior of 
the collective mode changes dramatically as the number 
of layers in the stack increases. Roton excitation become 
softened at a small but finite in-plane momentum. By 
contrast, when the magnetic field is parallel to the pan- 
cakes and q is along the field (i.e. = 7r/2 and a = 0, Fig. 
[TJb)), the energy of excitations is real and goes up with 
the increase of the number of layers. Dependence on the 
number of layers is not present for in-plane magnetic field 
but q being perpendicular to the direction of the field, 
since x(7r/2,7r/2) = in Eq. ([2]). These features can be 
understood using the following simple argument. Dipolar 
interaction is attractive when dipoles are oriented head- 
to-tail and repulsive when they are oriented side-by-side. 
So for perpendicular magnetic field increasing the num- 
ber of layers increases the attractive inter-layer compo- 
nent of the interactions but does not affect the repulsive 
intra-layer part. Conversely for the field parallel to the 
plane of the pancakes, increasing the number of layers 
primarily increases the repulsive part of the interactions. 
In Fig. Hfd), we show the collective mode dispersion for 
various values of with ten layers, showing a roton soft- 
ening when as < — 0.471ao (ao is Bohr radius). One pos- 
sible method of observing roton softening in experiments 
is to use Time-of-Flight experiments to measure the oc- 
cupation of the in-plane momentum states. When is 
close to the critical value of roton softening, the ther- 
mally excited quasi-particles begin to occupy the roton 
minimum and an incoherent peak at roton wavevector 
emerges. After releasing atoms from the trap, atoms pro- 
pogate essentially as free particles thus quasi-particles at 
roton wavevector expand in the in plane direction much 
faster than the condensate particles. As a result, for a 
sufficiently long expansion time, these particles should 
form a ring structure in the x — y plane. In the insert of 
Fig. md) we show thermal particle occupation number, 
^th(q, ^) = [exp(cj(q, /c; 0)//cbT) — 1]~ , for L = 10 but 
different 5-wave scattering length. 

We now proceed to discuss the structure of unstable 
modes as a function of their momentum q for different 
values of the condensate momentum (ko). Typical re- 
sults for the perpendicular magnetic field are shown in 
Fig. [21 Several interesting features should be noted: (1) 
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FIG. 2: Imaginary part of the collective mode as a function 
of \q\d and kd for magnetic field perpendicular to the layer 
plane. The upper, middle, and lower panels are for as — 
1.88ao, — 0.48ao, and = — l.SSao respectively, kod = 0.47r, 
and O.Gtt respectively for the left and right columns. All other 
parameters are the same as in Fig. [TJd). 



When the 5-wave scattering length is positive, there is no 
instability for k^d < 7r/2 (see (a)). Dynamical instability 
occurs when k^d > 7r/2 (see (b)). (2) When the 5-wave 
scattering length is negative but small, there is an insta- 
bility due to the roton softening for kod < tt / 2 (see (c)). 
(3) When the 5-wave scattering length is negative and 
large, there is another instability channel coming from 
the usual 3D collapse (see (e)). (4) Finally, the compli- 
cated interplay between the dynamical instability, roton- 
softening, and the collapse takes place for kod > 7r/2, as 
can be observed in (d) and (f). Such abundant struc- 
ture of unstable modes is unique to dipolar condensates 
in multi-layer systems. We note that earlier theoretical 
studies of the dynamical instability [12] did not consider 
the in-plane character of excitations. Hence all the earlier 
analysis was performed for systems which were effectively 
one dimensional. Our results indicate that the in-plane 
structure of collective modes is crucial when analyzing 
systems with dipolar interactions. 

In Ref. [9], the decoherence rate of BO was de- 
fined based on the growth rate of the momentum dis- 
tribution of atoms in the z-direction. Our next step is 
to connect our analysis of instabilities of states with a 
fixed momentum ko to the dynamics of BO. This can 
be done using following reasonable approximations: (1) 
In the process of BO, ko undergoes periodic oscillations 



and covers the entire Brillouin zone. Since the num- 
ber of non-condensate particles grows exponentially with 
a rate proportional to the imaginary part of the col- 
lective mode [13], we can calculate the growth rate of 
the momentum distribution (7(q, A:)) in a BO by aver- 
aging Im[c(;(q, /c; /co)] over all possible values of ko: i.e., 
7(q, /c) = ^ Jo"" dko |Im[cc;(q, /c; /cq)] | . (2) In the experi- 
ment, condensates are prepared at a fixed positive value 
of the scattering length (ain) before the magnetic field 
is changed abruptly to a different value for studying BO 
[14]. In our calculation, we use Bose-Einstein distribu- 
tion for the initial non-condensed particle distribution, 
i.e. nin(q, k) = nth(q, k) with ain = Sao and T = 100 
nK. (3) To include the effect of in-plane confinement, 
we introduce an infrared in-plane momentum cut-off, 
which should be of the order of the inverse of the system 
size. In other words, we assume that excitations with 
|q| < Qc are not relevant. (4) Finally, to simplify cal- 
culations, we assume that the condensate density is the 
same in all layers . We believe that these approximations 
correctly capture the phenomena taking place in experi- 
ments reported in Refs. [9]. They also provide a general 
framework for understanding the role of interactions in 
BO experiments with dipolar condensates. 

The time dependence of the width of the momen- 
tum distribution can now be calculated: Kz{t) = 

TtiEj" J^riU^,k)Pe^^(^''')f'^\ where /'dq 
is the integral of all the in-plane momentum with |q| > 
Qc J^in = X ^fe / ' (2^)^ ^in(q-fe) IS the total number of 
initial non-condensate atoms, in the limit of small time, 
K,{t) = K,{0){l^n^---) with 



dq 

(2^ 



■nin(q,/c)fc27(q,fc) (4) 



being the growth rate of the cloud size in z direction 
(i.e. the decoherence rate defined in Ref. [9]). Eq. (??) 
establishes a connection between instabilities of collective 
modes and the decoherence rate of BO. 

In Fig. [3l we show the calculated F as a func- 
tion of ttg, for ten layers of ^^K atoms. Here we use 
qc = 1.33 /d ~ 27r/i?||, where R\\ = 2.42 /im is the oscilla- 
tor length (i.e. approximate system size for small a^) of 
the in-plane trapping potential with trapping frequency 
44 Hz [9]. We point out several important features: (1) 
for the magnetic field perpendicular (parallel) to the layer 
plane, F has a minimum at = — 0.5ao(+0.38ao). These 
are different from zero due to dipolar interaction effects 
as originally suggested in Ref. [9|. (2) The minimum F 
we obtained can be considerably smaller than 1 Hz. This 
suggests that the maximal coherence time observed in ex- 
periments may be limited by other mechanisms, such as 
inhomogeneous density between different layers or laser 
stability [9]. Such additional effects not included in our 
analysis may also explain why our calculated positions 
of minimum decoherence rate are of the same order but 
still different from the experimental results [9I]. (3) We 
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FIG. 3: Decoherence rate of a condensate BO as a function 
of as. The solid (dashed) hne is for = 0(7r/2). qcd = 1.33 
and other parameters are the same as used in Fig. [TJd). In- 
set: Magnified results near the minimum value. The red (or 
higher) curves are results using qcd = 0.8 for comparison. 

emphasize, however, that the qualitative features of the 
decoherence rate are in good agreement with experiments 
in a wide range of scattering length. For example, the 
decoherence rate is not symmetric for the positive and 
negative scattering length regime when away from the 
position of minimum F. This is not surprising since in 
the former case the decoherence is dominated by the dy- 
namical instability (for \ko\d > 7r/2) whereas for the later 
case it comes from the global collapse (for \ko\d < 7r/2 ). 
Furthermore, the dependence of the decoherence rate on 
the scattering length has a positive curvature near the 
minimum, and a negative curvature for larger positive 
values of . The magnitude of F in the large regime 
is of the order of tens of Hz, also in agreement with the 
experimental results of Ref. \9]. (4) Finally, we find that 
the decoherence rate is highly dependent on the initial 
momentum distribution, system size, and other system 
parameters. For example, in the inset we also present re- 
sults calculated with a smaller infrared cutoff {qcd = 0.8), 
i.e. a larger in-plane system size. We observe a consider- 
able increase in the value of decoherence rate (red lines). 



indicating the important contribution from the instabil- 
ity in the long- wavelength limit. Detailed analysis (not 
presented here) also indicates that the decoherence rate 
can be very sensitive to the initial system temperature 
as well as the 2D density of atoms within individual pan- 
cakes. All of these effects originate from the in-plane in- 
stabilities in the stack of multi-layer dipolar condensate, 
which were not included in earlier theoretical analysis. 
Therefore, our model provides an alternative mechanism 
for explaining the decoherence rate of BO oscillations in 
experiments by Fattori et al. We point out that our 
analysis provides a fundamental limit on the decoherence 
rate of Bloch oscillations which arises from dipolar inter- 
actions for ultracold atoms and molecules. 

Before concluding this paper we would like to men- 
tion that our results imply greater stability of supersolid 
phases in multi-layer systems. It was argued before that 
roton softening can lead to supersolid phases which cor- 
respond to macroscopic occupation both at zero and ro- 
ton softening wave vector. In the case of a single layer 
system, roton softening occurs at short wavelengths (or- 
der of layer width) , and therefore is expected to lead to a 
global collapse after exciting higher transverse modes. In 
a multi-layer system we consider here, the roton softening 
occurs at wavelengths much larger than the layer width. 
This strongly suggests that supersolid phases should be 
considerably more stable in multi-layer systems. 

In conclusion we analyzed theoretically collective ex- 
citations in a multi-layer stack of two dimensional dipo- 
lar condensates. We found strong enhancement of roton 
softening and discussed its interplay with the dynami- 
cal instability. We showed important consequences of 
mode softening for the decoherence rate of Bloch oscil- 
lations. Our results are in qualitative agreement with 
experiments in Ref. [9]. We also make several concrete 
predictions for future experiments. 
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